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Abstract. The magnetospheric convection electric field contributes to Birkeland 

currents. The effects of the field are to polarize the plasma by displacing the bounce 

paths of the ions from those of electrons, to redistribute the pressure so that it is 

not constant along magnetic field lines, and to enhance the pressure gradient by the 

gradient of the bulk speed. Changes in the polarization charge during the convection 

of the plasma are neutralized by electrons in the form of field-aligned currents that 

close through the ionosphere. The pressure drives field-aligned currents through 

its gradient in the same manner as in quasi-static plasma, but with modifications 

that are important if the bulk speed is of the order of the ion thermal speed; 

the variations in the pressure along field lines are maintained by a weak parallel 

potential drop. These effects are described in terms of the field-aligned currents in 

steady state, isotropic, MHD plasma. Solutions are developed by taking the MHD 

limit of two-fluid solutions and illustrated in the special case of Maxwellian plasma 

for which the temperature is constant along magnetic field lines. The expression ^ 

for the Birkeland current density is a generalization of Vasyliunas’ expression for 

the field-aligned current density in quasi-static plasma and provides a unifying «£ 

expression when both pressure gradients and ion inertia operate simultaneously as r- w 

sources of field-aligned currents. It contains a full account of different aspects of the Y ™ 

ion flow (parallel and perpendicular velocity and vorticity) that contribute to the o o 

currents. Contributions of ion inertia to field-aligned currents will occur in regions z 5 

of strong velocity shear, electric field reversal, or large gradients in the parallel 

velocity or number density, and may be important in the low-latitude boundary 

layer, plasma sheet boundary layer, and the inner edge region of the plasma sheet. 


Introduction 

Birkeland currents in the Earth’s magnetosphere con- 
stitute the mechanism by which the magnetosphere is 
coupled to the ionosphere. The drifts of charged parti- 
cles, electrons and ions, in the magnetosphere carry cur- 
rent perpendicular to the magnetic field lines. The per- 
pendicular currents generally cannot close in the mag- 
netosphere. The requirement that the currents must 
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close on an Alfv6n transit timescale in order to main- 
tain quasi-neutrality means that they must close in the 
ionosphere via Birkeland currents. The closure currents 
determine the electric field in the resistive ionosphere, 
which in turn determines the electric field in the mag- 
netosphere, modifying the magnetospheric drifts. The 
behavior and evolution of the system is thus determined 
by the Birkeland currents. The convection electric field 
is always present and plays a primary role in the deter- 
mination of the currents. Any description of the origin 
of Birkeland currents should in principle take the con- 
vection electric field into account. 

In the context of fluid or MHD equations the effects of 
the convection electric field are described in terms of ion 
inertia, the time rate of change of the ion mean velocity, 
or equivalently, in terms of ion vorticity, velocity shear, 
or velocity gradients. It is known [Hasegawa and Sato , 
1979; Ogino , 1986; Southwood and Kivelson, 1991; Bim 
and Hesse , 1991] that ion vorticity contributes to field- 
aligned currents. The importance of the effects of ion 
inertia in magnetospheric plasma can be surmised from 
flow of plasma within different regions of the magneto- 
sphere. The association of field-aligned currents [Coley, 
1983], auroras [Chiu and Gomey , 1983], and the bound- 
ary plasma sheet [Lyons, 1981] with convection reversals 
suggest an intrinsic relation between field-aligned cur- 
rents and magnetospheric flows. In the plasma-sheet 
boundary layer, flow velocities of the order of hundreds 
of kilometers per second in a layer about an Earth ra- 
dius thick are known to persist for hours [Hones et al, 
1976; Eastman et al, 1984, 1985]. The low-latitude 
boundary layer [Eastman et al, 1976, Eastman and 
Hones, 1979] shows velocity gradients that result from 
the braking of the magnetosheath flow, of the order of 
1000 km s" 1 per Earth radius, and can be expected 
to contribute to the field-aligned currents as a result 
of the gradients in addition to viscous transport [Son- 
nerup, 1980; Phan et al, 1989]. Flows near the inner 
edge of the plasma sheet, which must turn sharply near 
the shielding layer, can also be expected to show con- 
tributions resulting from the ion acceleration. McFad- 
den et al [1986] observe that convection of cold plasma 
above the field-aligned acceleration region can account 
for their observations of precipitation at the edge of an 
arc. Galperin and Volosevich [1989] have shown that 
the polarization drift affects field- aligned currents if the 
electric field varies on the scale of the ion gyroradius 
and argue that it plays a role in the formation of nar- 
row, elongated, auroral arcs. 

It has long been argued that the theory of field- 
aligned currents is inadequate because existing models 
treat only one generator (pressure gradients, ion iner- 
tia, or viscosity) at a time [Akasofu, 1984] (see also Sis- 
coe and Fedder [1993]). For example, the theories of 
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Vasyliunas [1970], Rostoker and Bostrom [1976], and 
Sonnerup [1980] are restricted to pressure gradients, 
ion inertia, and viscosity, respectively, as the source of 
the field-aligned currents. The development of a the- 
ory that simultaneously treats all three is an important 
matter for geospace global circulation models that at- 
tempt to combine the different regions of the magneto- 
sphere (magnetosheath, magnetopause boundary layer, 
magnetosphere, and ionosphere) into a single, coupled, 
unified model. The purpose of this paper is to develop 
the theory of one level of unification of the generators 
- the simultaneous treatment of pressure gradients and 
ion inertia.- Specifically, we propose to quantify the re- 
lations among the field-aligned current density at the 
ionosphere, that is, the Birkeland current density, and 
the convection electric field and pressure gradients in 
steady state, isotropic, MHD plasma. Starting with a 
two-fluid description of the plasma, we derive an ex- 
pression for the current density that explicitly accounts 
for convection. The expression is a straightforward gen- 
eralization of those of Vasyliunas [1970] and Tverskoy 
[1982] to include the effects of electric and gravitational 
fields. The expression for the field-aligned current den- 
sity assumes that equilibrium is established along the 
magnetic field and is therefore valid on a timescale long 
compared to the particle bounce time and the Alfven 
transit time, the time it takes for an Alfven wave to 
travel the length of a field line. The results can be ap- 
plied to long-lasting Birkeland currents, but not rapidly 
varying systems such as substorm currents. Both re- 
strictions can be removed, but only at the expense of 
algebraic complications. 

The phenomenological basis for the assumption of 
isotropic pressure is the observation that magnetospher- 
ic plasma is often isotropic to a good approximation 
[e.g., Stiles et a/., 1978]. At the same time, the plasma 
is also collisionless, at least in the sense of Coulomb 
collisions. We note that purely collisionless drifting 
plasma generally cannot remain isotropic because the 
drifts energize the parallel and perpendicular motions 
of the particles at different rates (see Hines , 1963). The 
essential assumption of our analysis is that the bounce 
motions can be described by the behavior of collision- 
less isotropic plasma but that isotropy is maintained 
by some mechanism, not explicitly treated, on a drift 
timescale; such mechanisms include fluid (firehose and 
mirror) instabilities [Notzel et al , 1985], chaos in the 
plasma sheet [Buchner and Zelenyi, 1989], and strong 
pitch angle scattering [Kennel and Petschek, 1966]. The 
manner in which we take this assumption into account 
is to assume that the equatorial distribution functions 
are isotropic but that the bounce motions are collision- 
less. The assumption of isotropy can easily be replaced 
by one of gyrotropy using the methods of Heinemann 
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and Pontius [1991] (hereafter called HP91). To include 
gyrotropic pressure in the present paper would merely 
serve to complicate the presentation of what turns out 
to be a relatively simple physical picture. The mathe- 
matical details of time dependence and gyrotropic pres- 
sure will be published elsewhere. 

The principal formal result of the paper is the gen- 
eralization of Vasyliunas 1 [1970] formula for the field- 
aligned current density in quasistatic plasma to include 
the effects of ion inertia, gravity, and parallel potential 
drops; the reader is referred to equation (69) and those 
that follow. The expression, which is based on the as- 
sumption of isotropic Maxwellian distribution functions 
in steady state, is somewhat more complicated than 
that for quasi-static plasma but should be useful for the 
computation of Birkeland currents. The principal phys- 
ical results of the paper are that the convection electric 
field has three main effects on Birkeland currents: it po- 
larizes the plasma by displacing the bounce motions of 
the ions from those of electrons, it adjusts the pressure 
to satisfy equilibrium along the magnetic field, and it 
enhances the pressure gradient by the gradient of the 
bulk energy. The polarization causes a charge separa- 
tion of ions and electrons that convects at the E x B 
drift velocity. Neutralization of the excess charge by 
electrons is accomplished by currents that close through 
the ionosphere, that is, Birkeland currents. The paral- 
lel pressure equilibrium is maintained by a weak paral- 
lel electric field that survives in MHD approximation. 
The pressure gradients contribute to the Birkeland cur- 
rents through their gradients in the same manner as 
in quasi-static plasma, but with modifications that are 
important if the ion velocity is of the order of the ion 
thermal speed. In addition, the gradient of the square 
of the ion velocity contributes, essentially by adding to 
the pressure gradient. 

The paper is divided into three sections. In the first 
section we develop expressions for field-aligned currents 
in the terms of a two-fluid description of the plasma, as- 
suming that the ion and electron pressures are isotropic 
and the system is in steady state. In the second section 
we take the MHD limit of the two-fluid solutions and 
illustrate the results in terms of isotropic Maxwellian 
distribution functions for which the temperature along 
the magnetic field is constant. We conclude with a dis- 
cussion of the main results. 

Two-Fluid Analysis 

Prom the perspective of fluid behavior of magneto- 
spheric plasma the origin of the field-aligned currents 
is the plasma momentum. Gradient-curvature drifts, 
the ion inertial drift, and the gravitational drift all con- 
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tribute to the momentum density of the plasma, and 
must be balanced by the J x B force in steady state. 
This balance is expressed mathematically by the mo- 
mentum equation for quasi-neutral plasma 

+ dxvV + pV<f> = -JxB, (1) 

at c 

where p is the mass density, V is the bulk velocity, V 
is the kinetic tensor in the proper frame of the plasma, 
<t> is the gravitational potential, B is the magnetic field, 
J is the electric current density, and c is the speed of 
light; units are cgs. The tensor V is defined by 

Vi, = £ [PaSij + m a n a (V ai - Vi) (V aj - Vj )} , (2) 

a 

where a denotes the species of particle (electrons and 
ions), m a is the mass of the species, n a is its number 
density, P a is its scalar pressure, and V a is its mean 
velocity [see Rossi and Olbert , 1970, p. 290]. The con- 
vective derivative is d/dt = d/dt + V-V and reduces to 
V*V in steady state. 

Equation (1) determines the component of the cur- 
rent density perpendicular to the magnetic field in terms 
of the other quantities, that is, in terms of the perpen- 
dicular drifts. By itself, it is insufficient to determine 
the total current density. The parallel current density 
must be determined by the requirement that the total 
current density be solenoidal on an Alfv6n transit time 
scale, 

V-J = 0, (3) 

in order to prevent rapid charging of the plasma [see 
Vasyliunas , 1970]. 

The solution of (l)-(3) for the field-aligned current 
density has been approached in MHD approximation 
by a number of authors. The difficulty in obtaining 
solutions is the behavior of the inertial current 


, pc _ dV 
3p = 


( 4 ) 


which does not easily lend itself to mathematical analy- 
sis. Rather than work in terms of MHD plasma, a nat- 
ural way to include inertial effects is to follow Hasegawa 
and Sato [1979] by writing separate equations govern- 
ing mass and momentum of the electrons and ions. A 
benefit of working with separate fluid equations is that 
they lead to analogs of the one-fluid equations in which 
the equilibrium plasma is described by a set of field line 
constants and which can be solved in much the same 
way as in quasi-static plasma. The resulting two-fluid 
solution can be reduced by appropriate approximations 
to the MHD limit. 
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The laws governing the two-fluid motions are mass 
conservation, 

^ + V. n 0 V 0 = 0, (5) 

and the momentum equation, 

m a n a - a - ^° -f VP« +m 0 n 0 V<£ = n a q a E+ — V„xB, 

Uv C 

( 6 ) 


where E is the electric field. Equation (5) is valid for 
each species of particle (electrons and ions) as long as 
there are no sources or sinks of each species, for ex- 
ample, ionization, recombination, or charge exchange. 
Equation (6) is valid for plasma in which the momentum 
exchange between species can be neglected and in which 
the pressure can be regarded as isotropic. Anisotropic 
pressure and viscosity have been neglected; they can be 
taken into account but only at the expense of algebraic 
complications. The assumption of isotropic pressure ne- 
glects the so-called finite Larmor radius effects, which 
are important in collisionless plasma on the scale of 
an ion gyroradius [Macmahon, 1965; Stasiewicz , 1989, 
1991]. The convective derivative, d a /dt ~ d/dt + V 0 «V, 
depends on the species. 

To develop (6) in a form amenable to analytic treat- 
ment, we express the electric and magnetic fields in 
terms of potentials, E = — grad$> — (l/c)dA/dt and 
B = curl A, where $ is the electrostatic potential and 
A is the vector potential. Then the momentum equar- 
tion can be written 


dU a 


dt 


_ 

» n 


xr a 


= - — VPa - V (<7o$ + m a <j> + %m a v*) , 


( 7 ) 


where II a = m a V a + (q a /c)A is the “canonical momen- 
tum” associated with the fluid velocity, and the “gen- 
eralized vorticity” 

_ c „ „ ^ m a c .. 

r a — — V xn a = B H u) a (8) 

Qa Qa 

is the curl of the canonical momentum, and u; a = 
curlV a is the ordinary fluid vorticity; the units have 
been chosen such that T a has the same dimensions as 
B. 

We seek solutions in which equilibrium is established 
along the magnetic field, and therefore neglect the par- 
tial time derivative. In this case the velocity perpendic- 
ular to r o is 


V a± = 




r a x 


-VPa 


n a 


4- V (g a $ + fn a <j> -+- ^Tn a V£') 


( 9 ) 
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and the mass conservation equation (5) reduces to 

V-n 0 V 0 = 0. (10) 

It follows from (7) that 

IV [— VP„ + V (<?„$ + m Q ^+Im 0 V a 2 )l =0. (11) 
n a 

Equations (9)— (11) are two-fluid analogs of the quasi- 
static isotropic one-fluid equations = (c/B 2 )BxVP, 
divJ = 0, and B-VP = 0 that lead in quasi-static 
plasma to expressions for the parallel current density in 
terms of gradients of the flux-tube volume and pressure 
gradients [ Vasyliunas , 1970; Heinemann and Pontius , 
1990 (hereafter called HP90)]. Like B, T a is divergence- 
free; one may therefore associate flux tubes with T a 
analogous to magnetic flux tubes. By exploiting the 
analogies, one can include the effects of ion inertia. 

A solution of (11) can be obtained by writing the 
number density and pressure in the forms 



•/.(£; C,)cK (12) 

P.(X) = [{£- q a $ - m a <p - £m 0 v ; 2 ) 3/2 

3771a J 

■/«(£; Cj)d£, (13) 

where £ represents particle total energy and f a (£,Cj) 
is the particle distribution function expressed as a func- 
tion of particle total energy and arbitrary constants, Cj , 
that are used to specify the equatorial distribution func- 
tion. These are appropriate for plasma that is isotropic 
in the proper frame of each species (i.e., the frame in 
which the velocity of each species vanishes). One can 
verify by direct substitution that these solve (11) if it 
is assumed that r a -VCj = 0. 

The condition IVVCj = 0 is related to the polar- 
ization of the plasma by the convection electric field. 
Its physical significance can be sought in terms of a 
Vlasov description of the bounce motions of the parti- 
cles. From the guiding-center theory of plasma motions 
[e.g., Northrop , 1963] it can be shown that, if the Ex B 
drift velocity has a curl, then the guiding center velocity 
is given by 

v° = “if [ B + ( vx V ^)x] + • • • (14) 

where V# = cExB/B 2 represents the E x B drift and 
the ellipsis includes the E x B, gradient, curvature, 
and g x B drifts. The second term, which is equiva- 
lent to the polarization drift, means that particles are 
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displaced in their bounce motions by a distance that 
depends on the perpendicular vorticity. At an intuitive 
level this means that distribution function is constant 
along field lines defined by B H- (m 0 c/ < 7 0 )(curl Ve)j_ as 
are the constants characterizing the distribution func- 
tion. A tedious calculation (not given) shows that finite 
gyroradius corrections lead to F a -VC,* = 0[(m Q c/q Q ) 2 ] } 
where T a is based on the mean velocity of the species 
including the “magnetization” contribution to velocity, 
V< w) = ~{c/n a q a ) curl (P a B/B 2 ). At this order of an 
expansion in powers of m a c/q a this result does not dis- 
tinguish between the E x B-drift and the perpendicular 
component of the mean velocity of the species. It does, 
however, indicate the physical basis of the solution of 
(11) by (12) and (13) lies in the polarization displace- 
ments of the particles in their bounce motions. More- 
over, since it will ultimately be used only in the MHD 
limit, which also does not distinguish between V e and 
V a j_, it is entirely adequate for our purposes to assume 
that Fa-VCj = 0. 

As an example of (12) and (13), consider a convected 
Maxwellian particle distribution 


/o,e(v) 


n °.e -m 0 (V-V a,e) a /27o,e 

(27 rTo, e /m a ) 3 / 2 


(15) 


where n a , c is the equatorial number density, T a , c is the 
equatorial temperature in ergs, and V a?c is the equa- 
torial velocity. Because of the assumed isotropy in the 
moving frame, this can be expressed in terms of the 
total energy and equatorial constants: 


fa(S\C j ) = 


n { 




(27rT 0>e /m 0 )3/2 


£ - <7a$e - ~ ^Tn tt V% e 


(16) 

The number density and pressure that follow from in- 
tegration of (12) and (13) are 


n a = n a ^e~ qa( * a ~* a - e ^ Ta - e 

(17) 

p p p — 9o(^a — 'i'a,e)/Ta,e 

r a — * a,e c i 

(18) 

where represents the “effective potential” 

defined by 

<7o^a = <7a$ + rn a <f> + \m a V 2 . 

(19) 


(In writing (17) and (18), the effects of the loss cone 
have been neglected so that the pressure can be assumed 
to be isotropic all the way to the ionosphere. While 
this assumption is incorrect in application to magneto- 
spheric plasma, the magnitude of the error involved is 
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small, a few percent; see HP91 for an explicit evaluation 
in quasi-static plasma.) The number density and pres- 
sure are explicit functions of the equatorial constants 
n 0 ,fi, <£e, V^e, and T a>c . The meaning of (17) and 
(18) is that they are the number density and pressure 
at a point that is connected by field lines of T a to the 
equatorial crossing point with values n a>e , <j> C) V a >e , 

and T a ,e- 

Equations (12) and (13) solve the parallel equilibrium 
problem and thereby determine the three-dimensional 
distribution of n a and P a for given potentials and equa- 
torial distribution functions. The distributions of n Q 
and P a in turn determine the plasma drift velocities. 
Differentiation of (13) and use of (12) (or (18) and (17)) 
to identify n a shows that the following differential rela- 
tions hold: dP a /d$ = - q a n ai dP a /d<f> = -m a n 0 , and 
dP a /d\V% = — m a n 0 . These differential relations show 
that the gradient of the pressure, obtained by taking 
the gradient of (13), is 




ft p 

-n 0 V (g 0 $ + m a <t> + %m a V?) + ^ 

j J 


VCj. 


( 20 ) 

Because the parallel gradient of each Cj vanishes, that 
is, r.-VC, = 0 + 0[(m a c/ q a ) 2 ], (20) satisfies parallel 
equilibrium (11). It then follows from (9) that the per- 
pendicular velocity is given by 


V 0 _l = 



( 21 ) 


This expression for the perpendicular velocity, which 
is expressed in terms of quantities that are constant 
along r a , can be used to determine the field-aligned 
currents in the plasma. It is in a form analogous to the 
current density perpendicular to the magnetic field in 
MHD plasma and can be solved for the total velocity 
in the same manner. The solution is (see HP91 for 
methodology) 


9 an„V 0 (x) = c VV. j X VC,, (22) 

3 


where 



dP a d£ 
dC j r a 


(23) 


and the label T on the integral sign means that the 
integration is taken along T a from the equator to x. 
This solution is appropriate for a north-south sym- 
metric magnetosphere-ionosphere system such that the 
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Birkeland currents at the northern and southern ends of 
the flux tubes are identical. In general, there is an addi- 
tional force-free current aB running from ionosphere to 
ionosphere. While the results can be derived for more 
general configurations, including open field lines, we as- 
sume a symmetric closed field line situation for simplic- 
ity and simultaneously neglect any force- free current. 
The divergence of (22) is zero and its perpendicular 
component satisfies (9), showing that (22) is a solution 
of (5) and (6) in steady state, as required. 

The current density is given by the sum of (22) over 
species: 

( 24 ) 

a 3 

The Birkeland current density, that is, the vertical cur- 
rent density at the ionosphere, determined by methods 
identical to those of HP91, is 

Jr,i = -c^2 %^r 0 , e .v 0>e v„j X V a , c Cj, (25) 

a j 1 

where T a and r ori i are evaluated at the equator and 
ionosphere, respectively, the subscript r means the ra- 
dial component of a vector, V aie is a two-dimensional 
gradient evaluated at the equatorial crossing point of 
the field line of T a , and V a j is evaluated by integra- 
tion from the equator to the ionosphere. The subscript 
e denotes the equatorial crossing point of field line of 
the generalized vorticity and the subscript o will de- 
note the equatorial crossing point of a magnetic field 
line (see Figure 1). 

To demonstrate the physical content of (22)-(25), we 
give an example in terms of Maxwellian plasma. It 
follows from (17) and (18) that the derivatives of the 
pressure with respect to the constants are dP a /dP a , c = 
Pa/Pa,e = Pa/p a ,e, dP a /d$ e = q a Tl a , dP a /d<f> e = 
m 0 n 0 , dPa/d±V? e = m a n a , and dP a /dT a , e = n a q a (^a 
—'&a,e)/T a ,e. Here we have taken P ajC and T 0>c to be 
independent constants; the equatorial number density 
is not independent but is related to them by P a ^ e = 
n 0 , e r flt< .. The charge flux density (22) is 

q a n a V a = c(v f — ^xVP a , e 

\ Jr Pa,c 1 a 

+ v f ^xV( ?8 $ e + m.^ + |^) 

Jr A a 

^ f Qa n a{^ a ^a,e) 

Jr -Fa,c Fa 

(26) 

If parallel electric fields, gravity, and ion inertia are ne- 
glected, this reduces to 


xVT c 


-• 


Qa n aV i 


H? 


xVP tt 


10 



(27) 


where the label B on the integral sign means that the 
integration is taken along the magnetic field from the 
equator, the subscript o refers to the equatorial crossing 
point of the magnetic field line; we have assumed that 
the number density is constant along the magnetic field 
in the absence of a parallel electric field. The perpen- 
dicular component is the usual perpendicular velocity, 
V 0 x = cExB/B 2 + (c/n 0 g tt B 2 )BxVP a . The current 
density is given by the sum over species 

J=sc ( v X^ xVPo+v X^ xV$0 )’ (28) 

where P 0 is the sum of the electron and ion pressures at 
the equator; 77 is the charge density and can be assumed 
to vanish in quasi-neutral plasma. The first term in (28) 
is Vasyliunas 5 expression for the total current density 
due to pressure gradients in isotropic plasma. 

Equations (24) and (25) are formally exact, that is, 
(24) satisfies (1) and (3) exactly (recall, however, that 
r a-VCj* = 0[(m a c/q a ) 2 \). However, they are not very 
useful for analysis or modeling of field-aligned currents. 
As can be seen by comparison of (26) and (27), the ef- 
fects of a convection electric field are not only the addi- 
tion of explicit terms involving the velocity, but changes 
in the bounce paths as indicated by the differences in 
integration paths and equatorial values. The difference 
of the bounce paths, which reflects the effect of the vor- 
ticity on the bounce motions of the ions and contains 
information about the polarization of the plasma, is rel- 
atively small. The departure at the ionosphere from the 
field line occupied by a particle at its equatorial crossing 
point is of the order of the ratio of the equatorial Ex B 
drift velocity to the equatorial gyrofrequency, about 10 
km for hydrogen and a 100 km for oxygen, smaller than 
desirable for numerical simulations. On the other hand, 
the smallness of the difference in the bounce paths al- 
lows the development of accurate approximations, to 
which we now turn. 

MHD Limit of Two-Fluid Solutions 

Introduction 

The problem of the field-aligned current density is 
normally approached from the perspective of the MHD 
momentum equation 

p— 7 r + VP + pV<(> = ij xB. (29) 

at c 
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If inertia and gravity are neglected the solution for the 
Birkeland current density is [ Vasyliuna3 , 1970] 

Jr,i = /fxV.P. (30) 

Here B r i is the radial component of the magnetic field 
at the ionosphere, the subscript o means evaluation at 
the magnetospheric equator, and V Q is a two-dimensional 
gradient in the equatorial plane. 

The essential feature that allows the integration along 
field lines to give (30) is the constancy of the pressure 
along the magnetic field, B-VP = 0. The inclusion of 
inertia would be an easy matter if there were an analo- 
gous simple quantity that was constant along the mag- 
netic field; however, there appears to be no simple field 
line constant in that case. To see this, note that it 
follows from (29) in steady state that 

B-VP+pB-V (iy 2 + <j>) -pB-Vx (Vx V) = 0. (31) 

In steady state MHD approximation, with curl- free elec- 
tric field related to the bulk velocity by E = -V xB/c, 
the third term is — pV- [V x (VxB)], and therefore the 
parallel equilibrium equation is 

B*VP + pB-V (iy 2 - Vl + 4>) + pV- (BVj|V x ) = 0. 

(32) 

Not only is the pressure not constant along the magnetic 
field, but there appears to be no simple generalization 
of the constant pressure relation in terms of other con- 
stants. One can show that an equivalent form of the 
parallel equilibrium condition is 

^b-vp+b-v Ay 2 + <£+v-v J yj| dej =o. (33) 

if P — P(p) along B. The integral of the parallel veloc- 
ity is evaluated by integration along the magnetic field 
from the equator to arbitary points along the magnetic 
field line. The derivation of (33) from (32) is straight- 
forward but tedious. It is done by using the divergence 
properties of the last term of (32). The only possibly 
obscure step is to integrate the right-hand side over the 
surface of a flux tube, use Gauss 3 law to convert it to 
a surface integral, substitute the E x B drift, and per- 
form the surface integral by the methods of HP90. The 
result is 

J rf^rV- (PVj|V x ) = ~f (y Vg d^j d$, (34) 

where 3> is assumed constant along field lines, the con- 
tour integral is taken around any right-handed loop in 
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the plasma, and the integral over di is taken from the 
equatorial plane to the surface defining the ends of the 
flux tube. Applying Stokes 1 theorem to the right-hand 
side of (34) one finds 

J <PxV-(BV j,V x ) = ~cjvj Vj|<tfxV$-nda. 

(35) 

Combining this result, evaluated at the end of the flux 
tube with n = B/B , with the results of the more trans- 
parent integrals of the remaining terms of (32) gives 
(33). 

MHD Limit 

The direct solution of (29) for the parallel current 
density presents insurmountable difficulties. Rather 
than continue with that analysis, we solve the same 
problem by taking the MHD limit of the two-fluid so- 
lutions for the current density. The method is to ex- 
pand (24) in terms of the particle mass to charge ra- 
tios. By retaining only lowest-order terms consistent 
with the MHD approximation we obtain solutions for 
the MHD current density. To do this, we reexpress (24) 
in terms of integrals along the magnetic field lines and 
constants evaluated at the equatorial crossing points of 
magnetic field lines rather than field lines of the general- 
ized vorticity, by expanding each to order e a = m a c/q a . 
The approximation of Cj is based on the fact that it 
is constant along field lines of the generalized vorticity, 
Ta-VCj = O(e^), (see discussion after (13)). By ex- 
pressing r a in terms of the magnetic field and vorticity, 
one can write this as the differential equation 


^ = -^fw.-VC J + 0(4 ) , (36) 

where l represents arclength along the magnetic field. 
To order e a the constant Cj appearing on the right-hand 
side, which refers to the equatorial crossing point of the 
field lines of the generalized vorticity, can be replaced 
by which refers to the equatorial crossing point of 
the magnetic field line. Integrating along the magnetic 
field from the equator one obtains 


Cy = C j>e « Cj,o -^jff (Ua ' VCj ' o) % + <D ^ ’ 

(37) 


(where the label B on the integral sign means that the 
integral is taken along the magnetic field to arbitrary 
points in the plasma). The integration constant on the 
right-hand side is determined by the fact that in the 
limit that the evaluation point approaches the equator 
the field lines of the generalized vorticity and magnetic 
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field are identical. In a similar manner, one can show 
that 


- J r dCj r„ 

(0) _rr± c [ / vv (oA & 

°’ J q« J B ^ B 


+ 0(4), 


where 


/O) _ f 9P a d£ 

aJ ~ J B dCj B ■ 


(38) 


(39) 


Substituting these results into (24), one finds for the 
current density J(x) = J^(x) 4* J^(x), where 


J (0) (x) = wfjxvc ii0 (40) 

3 a 


J (1) (x) = 



E 


• [vvfjxv ^(o, a .VC ji0 )| 

+V / B ( u, “- VV S)f XV ^'° 


• (41) 


These hold at arbitrary points, x, in the plasma; so- 
lutions can be determined, in principle, by construct- 
ing the values of the integrals and field line constants 
as functions of the coordinates by integration from the 
equator to arbitrary points in the plasma and then tak- 
ing their local gradients. This is the general form of 
expansion of the current density and serves as the basis 
for evaluation of particular cases (i.e., specific distribu- 
tion functions). 

Equations (40) and (41) can be used to determine the 
Birkeland current density. When evaluated at the iono- 
sphere, each integral appearing on the right-hand side 
is integrated along magnetic field lines from equator to 
ionosphere and therefore can be regarded as a function 
of Euler potentials a and f3 that determine the mag- 
netic field according to B = VaxV/? [see Stem, 1970]. 
By expressing J in terms of a and /?, one finds for the 
Birkeland current density 


4? = -<= j §#b..£5> 


V^xVoCj.o 


(42) 


7 ( 1 ) _ J^r.iln rn a c 


(0) 


de 



Solutions for the Birkeland current density can be de- 
termined in much the same manner as the current den- 
sity in three dimensions, the difference being that the 
integrals must be evaluated by integration along the 
magnetic field from the equator to the ionosphere. 

Maxwellian Particle Distributions 

Equations (40) and (41) are expressed in completely 
abstract forms that depend on the constants character- 
izing the particle distribution functions. To use them 
for the computation of currents, one must state the 
distribution functions as functions of the total energy 
and the constants. Assume, for example, that each 
equatorial distribution function is a convected isotropic 
Maxwellian (16). The number density and pressure of 
each species are to be determined by integration of (12) 
and (13). To perform the integration, the limits of in- 
tegration must be considered. Two natural cases are 
potentials that increase or decrease monotonically from 
the equator. Of these, we treat only the former. In the 
latter case the particles see a potential maximum at the 
equator and are presumably lost in the ionosphere. The 
monotonically increasing potential is by no means the 
most general case. Nonmonotonic potentials are possi- 
ble and the reader should be aware of this restriction 
on our results. For monotonically increasing effective 
potential the lower limit of integration is the local (i.e., 
at the location £ along the field line) value of the effec- 
tive potential; the upper limit is S = oo. The number 
density and pressure that result from the integration 
are given by (17) and (18). Substitution for the partial 
derivatives with respect to the constants (see the text 
introducing (26)) into (40) leads to the current density 

j<o, - c ?( v /^i xvp - 

+ v / ^ x v (**• + + 

+ v J g ° na( r°~*° ,o) f : XVT.,.) 

+ <?(e + ) (44) 

(cf. (26)). The integrals are along magnetic field lines 
from the equator to arbitrary points in the plasma; the 
subscript B on the integral sign will be omitted here- 
after. At the same level of approximation the velocity 
that follows from (22) summed over species is 

pV = cV J V$o + O (e+) , (45) 

that is, the dominant contribution has been assumed 
to be from the electric field gradient, with the pressure, 
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gravitational potential, and velocity gradients neglected 
as being of one higher order in e + . For weak parallel 
potential drops the perpendicular component is the Ex 
B drift velocity, Vj_ = V^. The parallel component of 
the velocity generally does not vanish [see Schindler and 
Bim , 1987]. 

Solutions that satisfy the condition of quasi-neutrality 
are determined by the condition that the ion and elec- 
tron number densities be equal, that is, n = n + = n_, 
where ± denote ions and electrons, respectively. Ac- 
cording to (17), this condition is to be satisfied by the 
solution of 




-e(«' + -* +lC )/T +lC 


n_ |C e 




(46) 


(where e is the charge of the ions, assumed singly ion- 
ized) anywhere in the plasma as long the effective poten- 
tials of both ions and electrons increase monotonically 
with distance from the equator. According to (38), the 
relation between n± |C and n ±y0 = n Q is 

n ±,e = «o =F J (wVn„) ^ 4- 0(e±), (47) 


where uj = curl V is to be determined from (45). The 
condition of quasi-neutrality is then 

e(fr + - * + , e ) eQE- - *_,,) 

T+, e + T_, e 



= J (w-Vno) ^ + 0(4), (48) 


where m = m+ + m_ is the sum of the ion and electron 
masses. The potentials on the left-hand side can be 
referred to the equatorial crossing point of the magnetic 
field line by writing 

e(*+ - tf+.o) e($r_ - 0 ) 

7+,c T_, e 



, eOI'+.e - *+, 0 ) , e(«_, e - *-,o) 

TZ + TZ 
+ 0 ( 4 )- ( 49 ) 

In the last two terms on the right-hand side we retain 
only the electrostatic potential, regarding gravity and 
kinetic energy gradients as one higher order in e. Then, 
by (38), they can be expressed 


e(*±,« ~ *±,o) _ m±c 


T ±I , 


= T 


iS/ (l 


us* V 


\ _ , o \ 

,)-g+c?(4). 


(50) 
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By an analysis similar to that outlined in (34) and (35), 
one can show that 


J («.V*o)f = -~ c ( V 2 - v 0 2 ) + iv.v f Vj \dL (51) 


The solutions of (49)-(51) for the effective potentials 

are 


±eQE± -»±, e ) 

T±,e 


= ^[-1(^-1?) + (*-<« 

+ V.v/v l< «] + 0(« + ), (52) 


where T 0 is the mean value of the electron and ion tem- 
peratures at the equator. By (17) and (18) the number 
density and pressure are 


n = n 0 exp j ~ ^ ( V 2 - V?) 

'J V,|J}+0( e+ ) 


-v-v 


(53) 


P = Po exp{^- (V 2 - V 0 2 ) -(<£- *.) 

-V.vy , Vj | d«]|+0(e + ) l (54) 

valid if -V 2 /2 + <£> + V-V / Vj] c££ increases mononton- 
ically with increasing distance from the equator (see 
(52)). Differentiation of (54) along the magnetic field 
leads to (33), showing that (53) and (54) satisfy MHD 
parallel equilibrium (since, by (53) and (54), P = P(p)). 
The associated electrostatic potential drop is 


e($-S„) 
T-, o 


-V-V J V,di\ +0(e + ). (55) 


(The integral, as usual, is determined by integration 
from the equator to arbitrary points in the plasma). 
The entropy per unit mass, which is useful in writing 
the results, follows from the addition property of the 
entropy, ps = J2 a n am a s a , or 


S = So + — 

T 0 


4 (v 2 - v o) 

+ fo-&) + V-V J V\,dt 4- O (e+) , 


(56) 


where 




(57) 
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and the ions and electrons have each been character- 
ized by three degrees of freedom. Recognizing that the 
charge density is zero in quasi-neutral plasma, the cur- 
rent density (44) can now be written 

j<o)=e [ v /£ f xvp " 

+V J £^xv(^ + iv 0 2 ) 

+ v/ p(s-s 0 )^xVT 0 

+ 0(e+). (58) 

The current density can be treated in an analo- 
gous manner. We retain only those terms proportional 
to gradients of the electrostatic potential, regarding 
those involving thermal properties and gravitation as 
being of a higher order in e a . Then J^Ms 

jd) = -E m “ c2 [ v / ! ^ xV / 

♦*7 (-•*/« 

-bO(c + ). (59) 


Since neither the number density nor vorticity depends 
on the species in MHD approximation, (59) can be 
summed to yield 


J(1) , -C 2 V y ^xV J (wV*„)^ 


+ V 


J (a 


xV$„ 


+ 0(^), (60) 


where 

M(x) =J^ (61) 

is the mass per unit flux, a three-dimensional function 
of the coordinates determined by integration along mag- 
netic field lines from the equator to x. The first term 
can be written using (51). Using the fact that the vor- 
ticity has no divergence, one can rewrite the second 
term: 

J ( W -VM)| = J V-(VxVM)| 

= J V-(V fi xV|M + VnxVxM 
di 

+ V^xViM)-. (62) 

The first term of (62) can be expressed in terms of the 
electric field and the last can be integrated immediately. 
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The result is 

/ (u ,.VM)f = -J V- (^E X - ViieuxVxM) | 

+ ^Ex-VM, (63) 

where we have used VyM = p/B. In ordinary dielec- 
tric media, characterized by the plasma low-frequency 
dielectric constant \+Aivp<? fB 2 , the polarization charge 
density is 



In the present case the change in the number per unit 
flux on a given magnetic field line due to the presence 
of an electric field, obtained from (38), is 


6 


f nodi = _( f node _ f n a d£ \ 

1 B - U r a J B B ) 


-¥/(W 


n a {t')dt! 

B(i') 



That is, in comparison with the situation in the absence 
of an electric held, the particles are displaced from the 
field line, and have a different flux tube content on the 
given field line. The change in the charge per unit flux 

£/^ = «/(w-VAO§, (66) 

and represents the polarization charge per unit flux in 
the presence of an electric field. Consequently, by anal- 
ogy with ordinary dielectric media, we define the flux 
tube average of the polarization charge density by 

fnd* sc f(u,'VM)%. (67) 


With these identifications the current density is 

j ( 1) = -cv j v (v; 2 - v 2 + v-v f vj| dej 

- cV j ^ XV* 0 + O («+) . (68) 

The total current density is given by the sum of (58) 
and (68): 


J = c [ v /£f ><vp " 

+ v J ^xv(-iy 0 2 + <k 

■+■ v 2 — V * V J vj| dej 
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(69) 


+ V J p ( s -s 0 )^xVT 0 

-vj 


B 


+ 0(e+). 


As previously, solutions can be determined by con- 
structing, by integration from the equator to arbitrary 
points in the plasma, the integrals and equatorial con- 
stants as three-dimensional functions of the coordinates 
and taking the local gradients. The Birkeland current 
density, obtained by mapping the radial current density 
at the ionosphere to equatorial properties as in (42) and 
(43), is 


Jr,i 


XV 0 P 0 


fp.de 
J PoB 

J ^ xV °(-iK 2 + *o 

+v?-v 4 .v, J 


= -ci^Bo-lv, 
+ Vo 


+v e 


-v e / 


J - *o)j 

•qpde 


xVoTo 


B 


•XVo$o 


+ 0(«+). (70) 


where V* and V* are evaluated at the ionosphere. The 
nature of the solution is similar to that given by (69) 
except that the integrals are to be constructed by inte- 
gration to the ionosphere. If the convection velocity and 
gravity are neglected, (70) depends only on the pressure 
gradient and reduces to Vasyliunas’ [1970] formula (see 
(30)). The new information contained is the effects of 
the convection electric field. 


Verification of the Formula 


A necessary condition for (69) and (70) to be a cor- 
rect description of currents in MHD plasma is that (69) 
has the same perpendicular current density as that of 
the MHD momentum equation (29) and simultaneously 
satisfy current continuity (3). That (69) solves (3) is 
obvious from the fact that its divergence is zero. In or- 
der to help convince the reader that our result is indeed 
a solution for the currents in MHD plasma, we system- 
atically reduce (69) to (29) (including the effects of the 
gravitational field). To do this, we write (69) in the 
form 

j = c(v J £^XVP 0 + V J yxv (^ 0 + \v?) 

+ v J p(5-5 0 )|xVT 0 

/ d£ 

(w-VSo) ^ 
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(71) 


-cV J (wVM)^xVij + 0(e + ), 

that is, as the sum of (58) and (60). The perpendicular 
current density is 


Jx = ^Bx (JxB) 


= ^ B x 


Xvp 0 + pv(iy 0 2 + ^ 0 ) 

IPo 


= +p{s-s a ) VT 0 -pvj (w.V*o) ~ 

+ c( a/'V$ 0 ) VM — c(w*VM) V$ c 
-j-0(e + ), 


(72) 


where we have used the facts that B-VC 0 = 0 for any 
field line constant C Q and 


B-V 


f Fdt 

Jir =F ™ 

for any F. Now, it follows from (45) that 

pVx(VxV) = -c(wV$ 0 )VM + c(w.VM)V$ o 
+ 0(«+). (74) 

Combining (51), (72), and (74), one obtains 


+pV (v 2 - \V 2 + <t> 0 - v.v J Vj \<x) 


J± 


+ p(s - So ) VT 0 - P V x (V x V) j 
+ 0(e + ). (75) 

The gradient of the pressure that follows from (54) is 

VP = —VP 0 

P O 

+ pV i (y 2 - v 2 ) - (4> - *.) - v-v J V t de] 

4* p(s — s 0 )VT 0 4- 0(€+), (76) 


where we have used P Q = 2noT 0 , the factor of 2 follow- 
ing from the definition of the number density as equal 
to that of either ions or electrons and the definition of 
the temperature as the mean value of those of ions and 
electrons. The perpendicular current density is there- 
fore 


J x = ^BX 


vp + ^v^y 2 
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— pV x ( V x V) + pV<j> 




= -^Bx[VP + p(V.V)V + pV^] 

+ 0 («+). ( 77 ) 


that is, the perpendicular component (69) is just equal 
to the MHD current density in steady state. 


Discussion 


In this paper we have developed a description of 
Birkeland currents in isotropic MHD plasma to include 
the effects of ion inertia. The results are valid for 
timescales long compared to the Alfv6n transit time and 
the particle bounce time. Beginning with a description 
of the currents in isotropic plasma that is collisionless 
in its bounce motions, we passed to the MHD limit. 
The resulting expression, which is specific to isotropic 
Maxwellian distributions for both ions and electrons, 
accounts for the currents due to ion inertia as well as 
those due to pressure gradients in steady state. 

The forms of the expressions for the current den- 
sity can be thought of as generalizations of Vasyliu- 
nas * [1970] and Tverskoy’s [1982] expressions. They 
have been integrated along the magnetic field and de- 
pend on gradients of magnetospheric plasma in much 
the same manner. The difference in detail, that deriva- 
tives of the local velocity, not just equatorial quanti- 
ties, appear in (69), appears to be harmless. While this 
complicates the computation of currents near the equa- 
tor, Birkeland currents (70) can be computed in the 
same manner as with Vasyliunas’ formula. For a given 
magnetic field model, for example, Tsyganenko [1987] 
or R. Hilmer and G.-H. Voigt (A magnetospheric mag- 
netic field model driven by dynamically varying phys- 
ical input parameters, submitted to Journal of Geo- 
physical Research, 1994) the different contributions to 
(69), including the polarization current density, can be 
mapped to the ionosphere to provide models of field- 
aligned currents. The computation is straightforward 
if it is assumed that the equatorial distribution func- 
tion and electrostatic potential are known and that the 
parallel velocity is not important. To a first approxi- 
mation, regarding the field lines as equip otential, the 
E x B drift velocity is then determined along the mag- 
netic field, along with the potential drop and the mass 
density and pressure. 

The expressions for the current density are somewhat 
more complicated than the similar expressions for the 
current density in quasi-static plasma [ Vasyliunas, 1970; 
HP90]. There are four separate contributions to the 
field-aligned currents: the pressure gradient, the gradi- 
ent of the kinetic energy per particle, the gravitational 
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drift, and convection of polarization charge. The com- 
plications arise because equilibrium along the magnetic 
field is not characterized by constant pressure, as in 
quasi-static plasma, but depends on the electric field. 

The parallel component of the polarization current 
density (the last term of (69)) represents the diversion of 
polarization charge. If the polarization current density 
is defined as 

J (p) = — cV J x V$ 0 , (78) 

then the parallel current density is 

.v J^ = BVb . v J^ 

(79) 

and is equal to the time rate of change of the polariza- 
tion charge density. The origin of the Birkeland current 
density associated with the polarization is the neutral- 
ization of polarization charge by electrons. The convec- 
tion electric field has two simultaneous effects in this 
connection: it polarizes the plasma by displacing the 
bounce motions of ions from electrons and it convects 
both ions and electrons at the E x B drift velocity. 
Changes in the convection electric field or the density 
along the drift path change the polarization charge in a 
given flux tube. To maintain quasi-neutrality, electrons 
must move to neutralize the total charge. They can- 
not accomplish the neutralization by the E x B drift, 
which is in a direction perpendicular to the polarization 
displacement. Instead, they must neutralize the excess 
charge by closing through the ionosphere. 

The effects of polarization are normally discussed in 
terms of the parallel component of the ion vorticity, 
o/j| = B-u;/j5. In general, there are separate contribu- 
tions due to the parallel and perpendicular components 
of the vorticity. The contribution of polarization to the 
Birkeland current density written in terms of the par- 
allel and perpendicular components of the vorticity is 

J[ P) = — BVjj.V J (£1 + w r VM) (80) 

The contribution of the parallel vorticity is familiar, al- 
beit in slightly different form, from the work of Hasega - 
wa and Sato [1979]. Physically, it corresponds to the 
change in the volume occupied by ions in a given mag- 
netic flux tube due to the presence of an electric field. 
The second term, explicitly neglected in their analysis, 
describes the effects of the perpendicular vorticity and 
describes the effects of the displacement of bounce paths 
of the ions from those of electrons. 

The present results provide a direct way of calculating 
within a numerical model sources of current due to gra- 
dients of the velocity which can reasonably be expected 
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to make significant contributions in boundary layers and 
near the inner edge of the plasma sheet. The appli- 
cation of the resulting expressions may be useful for 
the analysis of Birkeland currents that originate in the 
low-latitude boundary layer, the plasma sheet boundary 
layer, and possibly the inner edge region of the plasma 
sheet. In the low-latitude boundary layer the solar wind 
velocity brakes from a few hundred kilometers per sec- 
ond tailward to a few kilometers per second earthward 
with flow speeds of the order of the ion thermal speed. 
Quite aside from the effects of viscosity, the velocity 
shears are likely to generate significant Birkeland cur- 
rents. In the plasma sheet boundary layer, parallel flows 
of the order of the ion thermal speed are observed on 
closed field lines and, according to Schindler and Bim 
[1987], have a large gradient near the separatrix be- 
tween open and close field lines. Reference to (70) and 
(62) suggests that this region is a source of field-aligned 
current. A further possible source is near the inner edge 
of the plasma sheet. Westward ion drift paths from the 
tail adjoin eastward drift paths that circulate the Earth 
[see Erickson et al . , 1991]. Across this layer the electric 
field reverses and the density need not be continuous, 
making it a candidate for Birkeland currents that orig- 
inate in the polarization of the plasma. 
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Figure 1. Schematic illustration of the relation be- 
tween ion bounce paths and the magnetic field. The 
perpendicular electric field points outward, as indicated. 
Ions are displaced from the magnetic field in the di- 
rection of the perpendicular electric field and follow a 
bounce trajectory indicated by the dashed line. Elec- 
trons, because of their small mass, essentially follow the 
magnetic field, indicated by the solid line. The equato- 
rial crossing point of the ion bounce path is indicated 
by e and that of the magnetic field by o. 
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